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ABSTRACT 

Fermi measurements of the high-latitude 7-ray background strongly constrain 
a decaying-dark-matter origin for the 1-100 GeV Galactic positron anomaly mea- 
sured with PAMELA. Inverse- Compton scattering of the microwave background 
by the emergent positrons produces a bump in the diffuse 100-200 MeV 7-ray 
background that would protrude from the observed background at these energies. 
The positrons are thus constrained to emerge from the decay process at a typical 
energy between ~100 GeV and ~250 GeV. By considering only 7-ray emission 
of the excess positrons and electrons, we derive a minimum diffuse 7-ray flux 
that, apart from the positron spectrum assumed, is independent of the actual 
decay modes. Any 7-rays produced directly by the dark-matter decay leads to 
an additional signal that make the observational limits more severe. A similar 
constraint on the energy of emergent positrons from annihilation in dark-matter 
substructures is argued to exist, according to recent estimates of enhancement in 
low-mass dark-matter substructures, and improved simulations of such substruc- 
ture will further sharpen this constraint. 



Subject headings: cosmic rays — dark matter — gamma rays: diffuse background 
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1. Introduction 

The positron excess in Galactic cosmic ray positrons between 1 a nd 100 GeV, despite a 
history of conflicting results, appears to be conflrmed by PAMELA (lAdriani et al.ll2009al ). 
There are, of course, several possible astrophysical explanations. Nearby sources (e.g. su- 



pernova or pulsars) of positrons, for exam ple, which suffer fewer 
positrons because they are younger (e.g. 



Eichler fc Maor 2005 



1989; 


Hooper et al. 


2009 


) or decay ( 


Eichler 


1989: 



■osses than typ ical Galactic 



Profumd 120081) . contribute 



suggested possible source. 

Annihilation, however, encounte rs a number of difficulties , as it requires a substantial 
boosting eit her from clumping (e.g. iKuhlen &: MalyshevI |2009[ ) or from a Sommerfeld en- 
hancement (jArkani-Hamed et al. 20091). and would lead to intense ph oton emission from the 



Galactic-Center region (e.g. lEichler &: Maori l2005l : IZhang et al.l 120091) or extr agalactic back- 



ground radiation from the superposition of haloes (jProfumo fc Jeltemall2009f ). that appears 
to exceed current observational limits. 

Here we investigate dark-matter decay as the source of the excess positrons observed with 
PAMELA. Our purpose is not a comprehensive theory of the dark-mat ter decay, in particu 



lar no t the viability of leptophilic decay, which is one of the challenges ( lArkani-Hamed et al. 



20091 ). given that standard hadronization scenarios are already excluded by the lack of an 
excess in the cosmic-ray antiproton data flAdriani et al.ll2009bl ). The decay source function 
scales with the dark-matter density, not the density squared as does annihilation, and there- 
fore the limits on decayi ng dark matter fr om observations of high-energy 7 rays from, e.g., 
dwarf galaxi es are weak dEpigLet al.ll2009l), b ut diff use emission may provide more stringent 



constraints ( llshiwata et al. 



20091: 



Chen et al.ll2009l ). We calculate the diffuse galactic 7-ray 



emission of the excess leptons, presuming they arise from dark-matter decay. The 7-ray 
intensity thus derived is more model-independent and also a lower limit to the true emission 
level, because we do not count 7 rays that are directly produced in the dark-matter decay, 
possibly via other unstable particles. 



2. Positron production by dark-matter decay 

2.1. The differential number density of positrons 

The propagatio n length of pos i trons and electron s at energies > 50 GeV is Lp ~ 
2^/trioss ^ 600 pc ( jPohl et al.l l2003l : iGrasso et al.l 120091 ) . considerably less than the scale- 
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length of the dark-matter halo in the Galaxy, i.e. the source distribution of positrons and 
electrons from dark-matter decay. Here, k is the energy-dependent diffusion coefficient and 
Tioss is the electron/positron energy-loss timescale. To first order, the spatial redistribution 
of those electrons by diffusive and convective transport can therefore be neglected. The 
differential number density of positrons has, in good approximation, the same spatial profile 
as the dark- matter density, pdm, because in decay processes the positron source function 
obeys Qf,+ oc Pdm- The differential number density of electrons/positrons in the steady state 
is then given by 

N{E,^) = -J— I duQ^yAu) (1) 
\E[,E)\ Je 

Above 10 GeV electron energy, the electron energy-loss rate, E{E), is dominated by syn- 
chrotron and inverse Compton losses. The soft-photon density varies somewhat between the 
Galactic plane and the halo, reaching a maximum of 1.5 eV/ cm ^ about 300 pc above th e mid- 



plane and falling to 0.8 eV/cm^ at 5 kpc above the midplane (jPorter fc Strong] |2005| ). Here 
we use an average value Up^ = 1.1 eV/cm^ to calculate the energy losses. The magnetic-field 
strength is not well known, and, using a canonical constant value for it, 8 /iG, we find for 
the energy-loss rate 



(2) 



where we ignore the weak modifications arising from the Klein-Nishina corrections to the 
inverse-Compton scattering rate. The parameter ( permits a scaling of the energy-loss rate 
relative to that given by the model for the galactic radiation field and a magnetic-field 
strength B = 8 /iG. 

2.7 eV /cm-^ 

The microwave background implies a lower limit ( > 0.1, and a magnetic-field strength 
B = 5 fiG would be described by C = 0.64. Estimates of the Galactic magnet ic field based 



on rn odeling the radio synchrotron emission give amplitudes from B = 6.1 fiG (IStrong et al. 



20001 ) to -B = 10 jjG ( jPohl fc Espositdll998l ). and the most likely value for ( is therefore in 
the range 0.75 < C ^ 1-3- The differential number density of electrons/positrons cannot have 
a spectrum harder than E""^, which is the spectrum that comes from the hardest possible 
injection spectrum. 

The PAMELA collaboration has measured a positron fraction that is rapidly increasing 
above 10 GeV. The highest-energy data point is 

+ 

'82 GeV) = 0.137 ±0.045 (4) 
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about 90% of which is in excess to the positrons expected from secondary production in 
cosmic-ray interactions. The total electron spectrum has not been determined with PAMELA 
yet, but we may use the excellent Fermi-LAT data th at between 20 GeV and 1 TeV are well 



represented by a single power law (lAbdo et al.ll2009l ) 



/ F \ -(3-045±0.008) 

Ne+ + Ne- = (73.5 ± 2.6) • 10"^ ( - — ) GeV'^ m^^ (5) 

\GeV / 

Combining Eqs. H] and [5] we find the positron density at 82 GeV as 

iVe+ (82 GeV) ~ 1.5 ■ 10^^^ GeV"^ m^^ (6) 

The positron excess over ordinary secondary production is measured to be only about 0.03 
at 17 GeV, thus indicating that the excess component has a spectrum not significantly softer 
than E'"^. Equation [T] indicates that this requires the injection of electron and positrons 
with a typical energy above 80 GeV. For simplicity we assume the injection of monoenergetic 
electrons, 

Qe+/e- = QoS (E — -Emax) (7) 

where -Emax ^ 500 GeV to avoid a spectral feature in the total electron spectrum, which 
would be in conflict with the power-law fit to the Fermi data. We shall see that the 7- 
ray limits calculated below provide much tighter constraints on -Emax- Assuming an equal 
number of excess electrons and positrons, we therefore estimate the total differential density 
of electrons/positrons that may come from dark-matter decay as 



E 



N{E, x) = (2.2 ■ 10-« GeV-i m"^) i?(x) J (E^^x - E) (8) 

In Figure [1] we demonstrate that the sum of this modeled dark-matter decay component 
and the expected contribution of secondary positron production provides a good fit to the 
PAMELA data. For comparison, the figure also displays the positron fraction for the case of 
a fiat injection spectrum, which may result from the decay of an intermediate particle with 
high kinetic energy. 

Qalt. = (^max ' E) (9) 



(-^max 


-E) 




' E 


/)( 


,GeV 



AT^it. {E, x) (X /?(x) \ \ - J J e (E^a. - E) (10) 

Here Qo is chosen 20% larger than in the case of monoenergetic injection to improve the fit 
of the positron fraction. To be noted from the figure is that -Emax ^ 200 GeV is needed for 
flat injection to well reproduce the positron fraction measured with PAMELA. 



0.01 



10 



E [GeV] 



100 



Fig. 1. — The measured positron fraction between 10 GeV and 82 GeV shown in compar- 
ison with a calculation for pure secondary production of positrons ( iMoskalenko &: Strong 



19981 ). the dark- matter contribution according to Eq. [8] (dotted line), and the total calcu- 



lated positron fraction (solid line). For comparison, we also show the positron fraction for a 
flat injection spectrum up to 200 GeV (dot-dashed line) and the total positron fraction for 
that case (triple-dot-dashed line). 
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i?(x) denotes the spatial distribution normalized to the local density, which should 
reflect that of the dark-matter particles. Simulations suggest that the dark-matter halo 
in the Galaxy is not spherical, and the vertical scale height may be only 40% of that in 
the Galactic plane fauhlen et al. 2009 ). The tidal streams pr oduced by the S agittarius 
dwarf speroidal may even indicate a triaxial density distribution ( iLaw et al.l 120091 ) . Here we 
conservatively assume a dark-matter halo with vertical flatteriing a nd for simplicity a global 



density scaling instead of, e.g., a NFW halo ( iNavarro et al.lll996l ). In that approximation, 



the normalized density proflle at the solar circle {r = Rq = 8 kpc) is 



i?(x) 



p{Ro,z = 0) 



1 + 2.5-^ 



2.2. The dark-matter lifetime 

Using Eqs. [H [21 and[H]we can estimate the source function of pairs from dark matter, and 
hence the dark-matter decay timescale. Assuming for simplicity decay into monoenergetic 
pairs, we flnd from Equations [1] and [7] 

Qo = C (5.5- 10~^° cm-=^s-^) (12) 

The source power is 

P = jdE EQ,y,- = QoE^,, = C (1.1 • 10^'' GeVcm-=^s-i) (^5^^) (13) 

Comparing with the local dark-matter mass density, po — 0.3 GeVcm"^, and accounting for 
an arbitrary branching ratio of conversion into pairs, rj, we flnd for the dark-matter decay 
time 

POV _.V ,r,r, , o26 ^ f ^max 



r..,^^.^(2.7.10^N)^^J (14) 

Care must be exercised in comparing astrophysical limits on the decay time on account of its 
dependence o n the energy of the injected pairs and on the assumed magnetic-fleld strength. 



For example, IChen et al.l ( 120091 ) use B = 3 pG, implying ( = 0.5, and consequently their 



limits on Tdecay are relatively high. 
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3. Inverse Compton scattering 



3.1. Galactic emission 



The electrons and positrons in the halo will inverse-Compton scatter soft photons into 
the 7-ray band where they can be observed with, e.g., the Fermi-LAT detector. The target 
photon field includes the microwave background, galactic infrared emission, and galactic 
optical emission, only the first of which is isotropic. The latter two will be backscattered 
toward the Galaxy, thus somewhat increasing the scattering rate com pared with the isotropi c 
case. The radiation field has been recently modeled out to z = 5 kpc (jPorter fc Strongj|2005l ) . 
Surprisingly, the energy density of optical light is higher in the halo than in the midplane 
on account of the thin-disk distribution of absorbers. 

The differential cr oss section for the up-scatt ering of isotropic photons of energy e to 
energy is given by (IBlumenthal &: Gouldlll970l ) 



a{E^,e,E) 



3 (Tt {rrie c 
4 



,2n2 r 



2glng + l + g- 2g^ + 



2{l + qT, 



where 



E^ 



and Fp 



AeE 
ml 



(15) 



(16) 



Fe {E — Ey 

The 7-ray intensity observed from the Galactic pole is then computed using the differen 
tial density of electron s/positrons (Eq. |8]) and the differential density of soft photons, n(e 
jPorter fc StrongI boosl ) . 



j^dz j dE N{E,x) j de n{e)a{E^,e,E) 



— I dz 

in 



(17) 



The resulting intensity is shown in Figure |2l To be noted from the figure is that even at 30 
GeV, where the emission arises from upscattering optical photons, the expected intensity is 
below 20% of the prelirninary estimate of the extragalactic isotropic background observed 
with Fermi ( lAckermanru 120091 ) , in line with the findings of llshiwata et al.l ( 120091 ) . Diffuse 
galactic emission from cosmic-ray interaction has an intensity similar to that of the extra- 
galactic isotropic background, and hence dark-matter decay is poorly constrained by 7-ray 
emission in the Galaxy. 

Whereas no significant intensity above 1 GeV is expected to come from beyond a distance 
of 5 kpc, the limit to which we have integrated the galactic emission, this is not true for 
the upscattering of microwave-background photons into the 100- Me V band. In fact, in the 
outer halo inverse-Compton scattering of the microwave background accounts for a larger 
share of the electron energy losses than near the Galactic plane, because the infrared and 
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optical photon fields quickly lose intensity beyond 5 kpc above the plane of the Galaxy. 
Also, the magnetic-field strength is expected to fall off, although we do not know at what 
point it drops to 3 /xG, below which the synchrotron energy losses are subdominant. We will 
estimate the intensity in the 100-MeV band from galactic dark-matter decay in the outer 
halo in the next section, together with the extragalactic component. 



3.2. Extragalactic background from intergalactic dark matter 

The majority of dark matter is located sufficiently far away from galaxies (> 20 kpc) 
that positrons and electrons from its decay would primarily interact with the microwave 
background. Since all electrons suffer the same fate and the electron source rate scales 
linearly with the dark-matter density, we can ignore any density structure and use spatially 
averaged quantities, i.e. a conservative fraction 90% of Qm = 0.239 times the critical density, 
9.9 • 10~^° gcm~^. For the source function and the energy-loss rate we find 

g,+/,- ~ (1 + zf C (2.2 • 10-=^^ cm-3 s-^) 5 {E - E^^) (18) 

and 



£;=(l + ^)V2.5-10-^^GeVs-^) (^^J 



(19) 



where z denotes the redshift. The average differential number density is therefore 



N{E, z) = ^ (8.8 • 10-^^ GeV-^ cm-=^) (^^^ O {E^^ - E) (20) 
The intensity of inverse-Compton emission is calculated as 

A = 4^(/^^^2 j^'^^f^^^ JdEN{E,z) J den{e)a[E,{l + z),e,E] (21) 

where 

f{z) = Jjj^ [(1 + {l + n^z)-z{2 + z) Qa] (22) 

is the Jacobian used to turn the line-of-sight integral into a redshift integral. We integrate to 
redshift z = 5, beyond which the cmissivity is negligibly small. The soft-photon distribution 
is the cosmological microwave background 
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with en = 2.4 • 10 ^ eV , added to which is a small contribution in the infrared band 



( iFranceschini et al.l Il998[ ) . The resulting expected 7-ray intensity is plotted in Figure [2] 
for C = 1) together with that produced in the outer parts of the Milky- Way halo at height 
z > 10 kpc where IC scattering off the CMB is also dominant and the expected differential 
density of electrons and positrons from dark-matter decay is 

N{E, x) = (2.4 ■ 10-^ GeV-i m-^) /?(x) (^^^ 6 (E^^x - E) , (24) 

larger than the galactic-plane solution (Eq. |H]) by the ratio of the energy-loss lifetime in 
the outer halo to that in the galactic plane. If the excess electrons and positrons could 
freely stream in the outer halo, we would be overestimating their 7-ray emission which 
contributes at most 1/3 of the total dark- matter signal. However, a number of well-known 
plasma instabilities would sharply impede such free streaming. As a 100-GeV electron has 
a propagation length of only ~15 kpc even if the diffusion coefficient is a 100 times that in 
the galactic plane, the effect of particle propagation is neglected in deriving Eq. [211 



At 200 MeV 7 -ray energy, the predicted intensity is close to that observed with Fermi 



(lAckermannl l2009l ). The uncertainty in the measured intensity is typically 15% and pre- 
dominantly systematic in origin. Because the electron/positron density is linear in (, the 
predicted 7-ray intensity would exceed the observational limits if the total magnetic-field in 
the solar vicinity were stronger than 10 /iG. It would also exceed the observational limits if 
the characteristic energy of the injected pairs were higher than ~250 GeV. For comparison. 
Figure [2] also shows the expected intensity for -Emax = 300 GeV to be twice that observed at 
200 MeV 7-ray energy. 

The same scaling with energy -Emax results for flat injection (cf. Eqs. [9l and [TO!) . The 
7-ray peak will appear at slightly lower energy and with somewhat reduced flux compared 
with monoenergetic injection at the same -Emax, because the mean particle energy is i?inax/2. 
This is demonstrated as well in Figure [2] where we show the expected 7-ray bump resulting 
from fiat injection up to 300 GeV, i.e. with characteristic energy ~ 150 GeV. 



4. Conclusions 



We have investigated 7-ray constraints on the notion that dark-matter decay is respon- 
sible for the recently measured excess of cosmic-ray positrons in the 1-100 GeV band. While 
the predicted GeV-band intensity from the inverse- Compton scattering of infrared and op- 
tical photons is below a prelimina ry estimate of the extragalactic 7-ray background based 
on Fermi data (lAckermannl |2009| ). the extragalactic background from the decay of inter- 
galactic dark matter would produce a bump at 100-300 MeV that is close to the observed 



0.01 0.10 1.00 10.00 100.00 

[GeV] 



Fig. 2. — Comparison of the predicted 7-ray intensity observed from the galactic poles with 
the preliminary estimate of the extragalacti c isotropic background observed with Fermi, 
indic ated by the das hed line (|Ackermannll2009r). and older data from COMPTEL ( IKappadath 



19991) and EGRET fISreekumar et al.lll998l : IStrong et al.ll2004l ). The thick line labeled Halo 



shows galactic emission according to equation [T71 and the other lines are for 7-ray emission 
from intergalactic dark matter, all for C = 1 and except where noted for E'max = 200 GeV. 
In detail, the thin solid line is extragalactic emission as given in equation [211 the dot-dashed 
line is emission from the outer Milky- Way halo beyond 2; = 10 kpc, the thick solid line the 
sum of the two, and the dotted line is the same for i?max = 300 GeV. The triple-dot-dashed 
line labeled FI is for flat injection up to 300 GeV (cf. Eqs. [9] and [TO|l . 
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extragalactic backgroun d at these energies, even one estimated from the older EGRET data 



( ISreekumar et al.l Il998[ ). Dark- matter decay therefore does not seem to be a viable ex- 
planation of the positron excess, unless the characteristic energy of the pairs produced in 
dark-matter decay is only in a narrow window between the lower limit ~ 100 GeV imposed 
by high-energy limit of the Pamela measurement and the high energy limit of ~ 250 GeV 
imposed by the extragalactic 100 to 200 MeV 7-ray background. 

Our stated limits are conservative in several ways: a) Any 7 rays directly produced in 
the decay of the dark-matter particles will lead to an additional signal that will make the 
observational limits more severe, b) We have neglected the contribution of astrophysical 
sources that are needed to explain the shape of the observed 7 ray spectrum, which is far 
broader than the dark-matter contribution. Any dark matter contribution would probably be 
noticeable at the 30 percent level, c) Any real decay process may involve multi-generational 
pairs, such as decay into electron positron pairs via taus and muons, even if quarks are 
avoided by some leptophilic process. This makes it harder to fit the pairs into the allowed 
energy window. 

We stress that the predicted 100-200 MeV 7-ray intensity from dark-matter decay is 
nearly model- independent, because it depends only on the total dark-matter density in re- 
gions outside those of strong galactic magnetic fields and starlight, i.e. in regions where 
inverse- Compton scattering of microwave background photons is a calorimeter for intergalac- 
tic electrons and positrons on account of its dominance among the energy-loss processes. The 
only assumptions we have made are time-independence of the dark-matter lifetime, which 
follows from conventional elementary particle physics, and the best estimates for the average 
magnetic field of the Galaxy in the solar neighborhood. 

Moreover, because the anomalous positrons reported by PAMELA are too energetic to 
retain their energy for a Hubble time, the only way to prevent the energy from going into 7 
rays is to preempt the inverse Compton losses with even faster synchrotron losses. It would 
seem then that magnetic fields of strength > 3 fiG are the only way to preempt IC losses. 

The argument would then apply to positrons from clump-enhanced dark matter an- 
nihilation as well as decay, if such enhancement occurred a) mostly in the outer part of 
galaxies where the magnetic field and starlight are weak, b) in extragalactic dark matter 
clumps or dwarf galaxies, where the magnetic field is weaker that in our own galaxy or c) 
elliptical galaxies, where the magnetic field strengths are likely to be less than in our own 
Galaxy. Quantifying this constraint for annihilating dark matter requires a reliable compu- 
tation of the distribution of annihilation, which is limited by numerical resolution at small 



mass scales. iDiemand et al.l (120071 ) argue that the positron emissivity per unit mass (q) in 



the outer halo, qoh, exceeds that in the local galactic neighborhood, qin, because substructure 
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survives tidal destruction far better there. This would raise the 7 ray contribution in the 
outer halo by a factor of qoh/lin relative to the limits we derive for decaying dark matter, 
where qoh/lin = 1 everywhere. Extragalactic substructure may give an even higher q, be- 
cause there is less tidal destruction, though there is the competing effect that in regions of 
lower density fluctuation, substructure is less likely to form at high redshift. Sommerfeld 
enhancement of the annihilation cross section at low center-of-mass energy would further 
raise the contribution of low mass substructures, in which the center-of-mass kinetic energy 
of the an nihilating pairs would be lower than in the diffuse dark rnatter in the solar neigh- 



borhood (IKuhlen. Madau. and Silk II2009I : iKistler fc Siegal-Gaskinsll2009l ). Thus, if we could 



make reliable simulations of dark-matter substructure down to the smallest mass scales for 
all cosmic locations, the results would probably already rule out explaining the PAMELA 
positron excess by annihilating dark matter. 

The arguments here do not, of course, immediately rule out detection of decay or anni- 
hilating dark matter by other more sensitive means, e.g. 7 rays from localized DM concen- 
trations. They merely argue against most scenarios for explaining the PAMELA positron 
anomaly in this manner. However, they call attention to the fact that any scenario for mak- 
ing 7 rays that includes accompanying e+ e- pairs, e.g. 7 rays by pion decay or final state 
radiation, may be constrained by the inverse Compton 7 rays that would accompany the di- 
rect 7 rays emitted in regions of weak magnetic field. Conversely, annihilation or decay that 
produces more pairs than direct 7 rays may be more readily detected via inverse Compton 
radiation of the pairs. 

We thank M. Kuhlen and P. Madau for helpful discussions. We acknowledge support 
from the U.S.-Israel Binational Science Foundation, the Israel Academy of Science, and the 
Robert and Joan Arnow Chair of Theoretical Astrophysics. This research was supported in 
part by the National Science Foundation under Grant No. PHY05-51164. 
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